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Abstract. We explore the connections between singular Weyl-Titchmarsh 
theory and the single and double commutation methods. In particular, we 
compute the singular Weyl function of the commuted operators in terms of 
the original operator. We apply the results to spherical Schrodinger operators 
(also known as Bessel operators). We also investigate the connections with the 
generalized Backlund-Darboux transformation. 



1. Introduction 

The present paper is concerned with spectral theory for one-dimensional Schro- 
dinger operators 

d 2 

(1.1) H = -— + q(x) 7 xe(a,b), 

on the Hilbert space L 2 (a,b) with a real- valued potential q E Ll (a,b). It has 
been shown recently by Gesztesy and Zinchenko [53], Fulton and Langer [13], |15j . 
Kurasov and Luger [33 that, for a large class of singularities at a, it is still possible 
to define a singular Weyl function at the basepoint a. Furthermore, in previous 
work we have shown that this singular Weyl function shares many properties with 
the classical Weyl function [30] and established the connection with super singular 
perturbations for the special case of spherical Schrodinger operators 

d 2 t 1(1 + 1) 
dx~ 2 

(also known as Bessel operators) |28j . 

On the other hand, commutation methods have played an important role in the 
theory of one-dimensional Schrodinger operators both as a method for inserting 
eigenvalues as well as for constructing solutions of the (modified) Korteweg-de Vries 
equation (see, e.g., [18] and the references therein). Historically, these methods 
of inserting eigenvalues go back to Jacobi [55] and Darboux [S] with decisive later 
contributions by Crum [7J , Krein [3Tj , Schmincke [5S] , and Dcift 9 . Two particular 
methods turned out to be of special importance: The single commutation method, 
also called the Crum-Darboux method [7J , [5] (actually going back at least to Jacobi 
[26) 1 and the double commutation method, to be found, e.g., in the seminal work 
of Gel'fand and Levitan [TB]. For recent extensions of these methods we refer to 



(1.2) H = -— + ^^+q(x), x€(0,oo), 
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Krein [31] was the first to realize the connection between inverse spectral prob- 
lems and Cram's results. Namely, in |31j . the connection between the spectral 
measures of the original and transformed operators was established and then ex- 
ploited to characterize the spectral measures of Bessel operators in the case I G N 
(see also [13]). This idea has been subsequently used by many authors: See, for 
instance, E] [13l [24] and references therein. 

The purpose of our present paper is to continue the work of Krein and estab- 
lish the connection between the singular Weyl functions (and hence between the 
spectral measures) of the original and transformed operators for both the single 
and double commutation method. In particular, we will obtain an independent 
proof for the fact that the singular Weyl function of perturbed Bessel operators 
is a generalized Nevanlinna function. In addition, we investigate the connections 
with the generalized Backlund-Darboux transformation (GBDT) for a particular 
example. This method is a generalization of the double commutation method which 
it contains as a special case (cf. Subsection 15. ip . 



2. Singular Weyl-Titchmash theory 

We begin by recalling a few facts from [30]. To set the stage, we will consider 
one-dimensional Schrddinger operators on L 2 (a,b) with — oo < a < b < oo of the 
form 



(2.1) T = - +q 

ax 

where the potential q is real-valued and satisfies 

(2.2) qeLj oc (a,b). 

We will use r to denote the formal differential expression and H to denote a corre- 
sponding self-adjoint operator given by r with separated boundary conditions at a 
and/or b. 

If a (resp. b) is finite and q is in addition integrable near a (resp. b), we will say 
a (resp. b) is a regular endpoint. We will say r, respectively H, is regular if both a 
and b are regular. 

We will choose a point c G (a, b) and also consider the operators HR C \, 
which arc obtained by restricting H to (a, c), (c, 6) with a Dirichlet boundary con- 
dition at c, respectively. The corresponding operators with a Neumann boundary 
condition will be denoted by n and HP . 

Moreover, let c{z,x), s{z,x) be the solutions of tu = zu corresponding to the 
initial conditions c(z,c) = 1, c'(z,c) = and s(z,c) = 0, s'(z,c) = 1. 

Define the Weyl m-functions (corresponding to the base point c) such that 

u-(z,x) — c(z,x) — m-(z)s(z,x), z e C\ a{HP ac ^), 

(2.3) u + (z, x) — c(z, x) + m + (z)s(z, x), z G C\a(HR Vj ), 

are square integrable on (a, c), (c, b) and satisfy the boundary condition of H at a, 
b (if any), respectively. The solutions u±(z, x) (as well as their multiples) are called 
Weyl solutions at a, b. For further background we refer to [121 Chap. 9] or [415] . 

To define an analogous singular Weyl m-function at the, in general singular, 
endpoint a we will first need an analog of the system of solutions c(z, x) and s(z, x). 
Hence our first goal is to find a system of real entire solutions 9(z, x) and 4>(z, x) 
such that 4>(z,x) lies in the domain of H near a and such that the Wronskian 
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W(9(z), (f>(z)) = 1. By a real entire function we mean an entire function which is 
real-valued on the real line. To this end we start with a hypothesis which will turn 
out necessary and sufficient for such a system of solutions to exist. 

Hypothesis 2.1. Suppose that the spectrum of HP c s is purely discrete for one 
(and hence for all) c G (a, b). 

Note that this hypothesis is for example satisfied if q(x) — > +oo as x — > a (cf. 
Problem 9.7 in [42]). 

Lemma 2.2 ([30]). Suppose Huvothesis \2. 1\ holds. Then there exists a fundamental 
system of solutions (f>(z,x) and 9(z 1 x) of tu = zu which are real entire with respect 
to z such that 

(2.4) W{6{z),<t>{z)) = l 

and(j)(z,.) is in the domain of H near a . Here W x (u, v) — u{x)v'(x) — u'(x)v(x) is 
the usual Wronski determinant. 

It is important to point out that such a system is not unique and any other such 
sytem is given by 

6(z, x) = e- 9{z) 8(z, x) - f{z)<f>{z, x), <j>(z, x) = e a ^(f>{z, x), 

where g(z), f[z) are real entire functions. 

Given a system of real entire solutions 4>{z, x) and 6{z, x) as in the above lemma 
we can define the singular Weyl function 

(2 .5) mm . -»^TT 

W{4>{z),u + (z)) 

such that the solution which is in the domain of H near b (cf. ([273])) is given by 

(2.6) u + (z,x) = a(z)(6(z,x) + M(z)<f>(z,x)), 

where a(z) = —W((j)(z),u + (z)) — f3(z) — m + (z)a{z). By construction we obtain 
that the singular Weyl function M(z) is analytic in C\R and satisfies M(z) = 
M(z*)*. Rather than u + (z,x) we will use 

(2.7) i>(z,x) = 6{z,x) + M(z)4>(z,x). 

Recall also from [30] Lem. 3.2] that associated with M(z) is a corresponding 
spectral measure 

1 1 f Xl 

(2.8) - (p{(x 0) xi)) +p([xo,Xi})) =lim- / Im(Af(x + ie))da:. 

A e\a 7T J Xg 

3. Connection with the single commutation method 

3.1. Preliminary basic results. We begin by recalling a few basic facts from the 
single commutation method. Let A be a densely defined closed operator and recall 
that H = A* A is a self-adjoint operator with ker(£f ) = ker(A). Similarly, H = AA* 
is a self-adjoint operator with ker(Jf) = ker(^4*). Then the key observation is the 
following well-known result (see e.g. [H] Thm. 8.6] for a short proof): 



Theorem 3.1 ([9]). Let A be a densely defined closed operator and introduce 
H = A* A, H 
ily equivalent. 



H = A* A. H = AA* . Then the operators H\) SBt (jg\± an d -^Ikerfff)- 1 are un ^ ar 



4 



A. KOSTENKO, A. SAKHNOVICH, AND G. TESCHL 



// Hip = Xip, X € R, tp e T){H), then $ = Aip e D(H) with Hip = Xip and 
= vAHV'll' Moreover, 

(3.1) R A (z) D - (AR H {z)A* - 1) , R H (z) D - (A*R 6 (z)A- l) , 

where Rh(z) = (H — z)^ 1 denotes the resolvent of an operator H . 

3.2. Application to Schrodinger operators. In order to find such a factoriza- 
tion for a given Schrodinger operator H one requires a positive solution of the 
underlying differential equation. Existence of such a solution is equivalent to semi- 
boundedness of H and we will thus make the following assumption: 

Hypothesis 3.2. In addition to Hypothesis \2.1\ assume that H is bounded from 
below and limit point at b. Let A € R be such that H — A > 0. 

In particular, let <p{z,x), 9(z,x) be a fundamental system of solutions as in 
Lemma [2.21 and recall (cf. [HJ Lem. 9.7]) that the Green's function of H is given 

by 

(3.2) G(g,x,y) = \^ Z,X ^ Z,y ^ X ~ Vl 

\cf)(z,y)ip(z,x), x>y. 

Moreover, if H — A > we must have <p(\,x) > (as well as ip(\,x) > 0) for 
x e (a, b) possibly after flipping signs by [19l Cor. 2.4]. 
Now consider the operator 

A f f d 4_££i£) 

(3.3) J)(^) = {/ € L 2 {a, b)\f G ACi oc (a, b), a<j> f € L 2 (a, b)}. 

Here we use A^ and for the operator and differential expression, respectively. 
It is straightforward to check (cf. [HI Problem 9.3]) that A$ is closed and that its 
adjoint is given by 

A *f_* f * _ d cp'{\,x) 

= U e i 2 (a, 6)1/ G ACocfo, 6), a;/ e i 2 (a, &), 

(3.4) lim /(x) 5 (x) = 0,V 9 €%)}. 

If we also have 8(X, x) > we can also define Ag by using 9(X, x) in place of 4>(X, x). 

Lemma 3.3. Assume Hypothesis \3.2\ holds. Then H — X = A*iA§. If in addition 
9(X, x) > and t is limit point at a, then we also have H — X = AgAg. 

Proof. It is simple algebra to check that H and A^A^ agree on functions with 
compact support in (a, b) hence A*^A^ is a self-adjoint extension of the minimal 
operator associated with r and it remains to identify the boundary conditions. 
Since by assumption r is limit point at b we only need to consider a. Moreover, 
since (p(X, x) is in the domain of A^A^ near a by our choice of A$, both A*^A^ and 
H are associated with the boundary condition generated by <f)(\, x) if r is limit 
circle at a. 

The case of Ag is even simpler since by our limit point assumption there is no 
need to identify any boundary conditions. □ 
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The commuted operator H — A = A^A*^ is associated with the potential 

Moreover, it is straightforward to check that if u(z) solves tu = zu then 

(3.6) u{z, x) = a,pu(z } x) = -— — r 

solves tu = zu and given two solutions u(z) and v(z) we have 

(3.7) W(u(z),v{z)) = {z- X)W(u(z), v(z)). 
Similarly, H — A = AoA* e is associated with the potential 

,„ „s / \ / n n d 8' (X, x) 

and 

(3.9) u(z, X ) = -a e u( Z , X ) = ^^m 

0(X, x) 

solves fu — zu and given two solutions u{z) and v(z) we have 

(3.10) W(u(z),v{z)) = (z- X)W(u(z), v(z)). 



Theorem 3.4. Assume Hypothesis \3.2\ holds. Then the operator H = A^A*^ + A 
has an entire system of solutions 

mil i, s 1 v W x (0(\)^(z)) f*<j>(X,y)<P(z,y)dy 

(3.11) J>(z,x) = j-jatfax) = - {z _ mXtX) = ^ > 

(3.12) e { z, x ) = a,e {z , x) = - w ^^\ 

which satisfy W(6(z),(f)(z)) = 1. In particular, H satisfies again Hvpothesis \2.1{ 



Furthermore, the Weyl solutions of H are given by 
(3.13) 

<j»(z,x), ij>(z,x) = asi>(z, X ) = - _ §( z x ) + M( z )cj>(z,x), 

(f>{X, X ) 

where 

(3.14) M(z) = {z- X)M(z) 

is the singular Weyl function of H corresponding to the above system of solutions. 
Moreover, the associated spectral measures are related via 

(3.15) dp{t) = (t- X)dp(t). 
Proof. The first part is easy to check using (|3.7[) . To see 

(3.16) W x ^(X),cf > (z)) = (z-X) ( ^(X,y)<f>(z,y)dy 



note that both sides have the same derivative and are both equal in the limit 
x — » a. 

For the second part we set A = and fix z 6 C\K without loss of generality. 
Moreover we will abbreviate <p(X,x) — 4>q( X ), 4>(z,x) = <fi(x), and ip(z, X ) = tp(x). 



6 



A. KOSTENKO, A. SAKHNOVICH, AND G. TESCHL 



Now let / be some absolutely continuous function with compact support in 
(a, b) and /' € L 2 . Then invoking (|3~T|) one computes using ([274]) . (|3TT|) . (pTT2|) . 
and integration by parts that 



G{z,x,y)A%f{y)dy- f{x) 



1>{x) / m{o.%m)dy +a 0(x) / i>{y){a%f){y)dy 

J a J x 

\ ((a^)(x)^(x)f(x) + J f(y)M)(y)dy) 

M)(x)(-i/j(x)f(x)+ f f{y){a^){y)dy) - f{x)\ 



+ 



=-[zi){x) I f(y)4>(y)dy + z<j>(x) I f(y)tp(y)dy 



+ijj(z)(f)(x)f(x) - z<f>(x)i/)(x)f(x) - /(x 

rX rb pb 

=$(x) I <j>(y)f(y)dy + (j>(x) / 4>(y)f(y)dy = / G{z,x,y)f{y)dy, 



which shows that 

(3.17) G(z,x,y) 



cj>(z,x)ii>(z,y), x<y, 
4>(z,y)i>(z,x), x>y, 



is the Green's function of H since the set of / under consideration is dense. Since 
the Green's function is unique, it follows that <f)(z,x) and ip(z 7 x) are the Weyl 
solutions of H. Finally, existence of the entire Weyl solution <f>(z, x) is equivalent 
to Hypothesis [O by Lemma 2.2 in [35]. □ 

Note 

(3.18) ^,) = -^^ ) 'f (A)) = r ^'^ i $(A, x) 1 



</>(A,x) </>(A,x) ' ' 0(A, x)' 

where the dot denotes differentiation with respect to z. 

Remark 3.5. A few remarks are in order: 

(i) If t is regular at a and 9(z,x), 4>{z,x) are chosen to satisfy the boundary 
conditions 

9(z, a) — cos(a), 9 (z, a) — — sin(a), 
4>{z, a) = sin(a), tj> (z, a) = cos(a), 

then t will be again regular if and only if sin(a) ^ 0. Moreover, in this 
case we have 

9{z, a) = sin(a) , 9'(z, a) = cos(a) {z — A — 1 — cot(a) 2 ) , 
(p(z,a) = 0, 4>'{z,a) = sin(a). 
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(ii) If(j)(X,x) is the principal positive solution near a (i.e., if H is the Friedrich's 
extension of t, see [34] . [22j ). that is, 

dx 



(3 ' 19) 'a 4>{\X? 

then f will be limit point at a. 
(iii) Note that H has no eigenvalue at X: o~ p (H) n {A} = 0. 

Recall that N%° is the subclass of N K consisting of generalized Nevanlinna func- 
tions with k negative squares and having no nonreal poles and the only generalized 
pole of nonpositive type at infinity. In particular, M £ iV£° admits an integral 
representation 

(3.20) MM . (i + I (jh - TT?) ITTFF + 

where Ao = inf cr(H), k < K, I < 2k + 1, 

(3.21) aj G R, and / (1 + t 2 )- fc_1 dp(i) < oo. 



Without loss of generality we can assume that the representation (|3.20p is irre- 
ducible, that is, either k = or L(l + t 2 )~ k dp(t) = oo. For further definitions and 
details we refer to, e.g., [32], [HI §2.2-2.3], [12 App. C]. 

Corollary 3.6. Assume Hvvothesis \S. 2\ holds. Assume also that the functions M(z) 
and M(z) are connected by (|3.14j) . Then M S for some k € No if and only if 
M € N?° with 



(3.22) 




lim aT oo §M G [0, oo) 
MM 



lim„ 



ytoo (ITpT G [-00,0). 



Proof. The proof is similar to the proof of Lemmas 5.1 and 5.2 from [2 7) . Assume 
that M £ for some k e No. Using an irreducible integral representation (|3.20l) 
we will show that 



M(z) = (z- A)M(z) 



' ^ dp(t) . J- 

, i / ^rr + (Z — An) > CLjZ J 

t-z 1 + t 2 J (l + t 2 ) k v 

' j=Q 



where Aq = inf a(H), lies in N~°. Using 



1 t \ , 2 . 1 t- A zA+1 



7 7T^ )=(z + 1 ) 



t-z 1 + t 2 J v 't-zl + t 2 l + t 2 

we obtain 

M(z) = (1 + z 2 ) k+1 / + °° ( — d P® + y a 

where dp given by (|3. 15|) is positive since A > Ao by Hypothesis 13.21 

The converse implication M £ N~° M E N%° can be established analogously. 
Finally, the connection between k and k is straightforward from the following 

characterization of -functions: Given a generalized Nevanlinna function M in 
, the corresponding k is given by the multiplicity of the generalized pole at oo 
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which is determined by the facts that the following limits exist and take values as 
indicated: 

~ W^ 1 ( ' °° ] ' y" W^ 1 1 ' °° ) - 

□ 

Similarly, we obtain 

Theorem 3.7. Assume Hypothesis VS.'zX r is limit point at a, and let 9(X,x) > 0. 
The operator H = A$Ag — A has an entire system of solutions 

3.23 <p{z, x) = -a e (f)(z, x) = — — r , 

b^A, X) 

which satisfy W(6(z), 4>{z)) = 1. In particular, H satisfies again Hvpothesis \2.1{ 

Furthermore, the Weyl solutions of H are given by 
(3.25) 

4>(z, x), ip(z, x) = -a e ip(z, x) = y^—TT^TT^T- = x ) + M(z)(f>(z, x), 

(z — X)t)(X, X) 

where 

(3.26) M(z) = (z- A)" 1 M(z) 

is the singular Weyl function of H . The associated spectral measures are related 
via 

(3.27) dp(t) = (f - X^dpit) - M{X)dQ(t - A), 

where 0(i) =0 for t < and 0(t) = 1 for t > is the usual step function. Here 
M(\) = lim £l0 M(X-e). 

Note 

no* ^ 1 fin \ W x {d{X),e\X)) 

Remark 3.8. Again a few remarks are in order: 

(i) If 8(X,x) is the principal positive solution near b, that is, 



(3.29) 



dx 



0(X,x) 2 

then f will be limit point at b. (Note that the principal positive solution 
near a is 4>(X, x) since we assumed the limit point case at a.) 

(ii) Note that H has an eigenvalue at X unless 9(X,x) is the principal positive 
solution near b. 

(iii) Note that factorizing H using 9(X,x) = <f)(X,x)~ l will transform H back 

into H . In particular, (j)(z, x) = 4>(z, x) and 6{z, x) = 8{z, x). 

(iv) Clearly this procedure can be iterated and we refer (e.g.) to Appendix A of 
|20) for the well-known formulas. 
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3.3. Example: The Coulomb Hamiltonian. We can apply the single commu- 
tation method to the Coulomb Hamiltonian 

d 2 | *(* + !) 7 

dx 2 x 2 x ' 



(3.30) Hi = -— + -±, IG No, 



by setting 

tt«m Af- d i 1 + 1 i 7 A*f- d i 1 + 1 i 7 

which gives (01 Thm. 10.10]) 

(3.32) Hi = A^Ai - cf, H l+1 = A X A\ - cf, 

where 

^ ci = wrvy 

In particular, the singular Weyl function is given by 

l-i 

(3.34) M l (z)=M (z)H(z-cl), 

k=0 

where the singular Weyl function M (z) of H is a Herglotz-Nevanlinna function 
since Ho is limit circle at a = 0. 

Remark 3.9. The singular Weyl function for this case was explicitly computed in 
p~5l eq. (5.11)]. Moreover, the fact that it can be reduced to the case I — via the 
above product formula was also first observed in [151 Lemma 5.1]. In the special case 
7 = we obtain Mi(z) = y/—zz l as was first observed by |llj (see also Section 5 in 

3.4. Application to perturbed Bessel operators. Next we want to apply these 
results to perturbed spherical Schrodinger equations which have attracted consid- 
erable interest in the past p], [5], [25] [29], [28]. In particular, we want to mention 
[2], where commutation techniques were used to transfer results for I — to I € Nq. 



Lemma 3.10. Fix I > —\ and p £ [l,oo]. Suppose 

d 2 t 1(1 + 1) 
dx 



(3.35) H = -—+ 2 +l( x )i ^(0,6), 



where 
(3.36) 



xq(x) G L p (0,c), p£ (l,oo],Z > -|, 

i^ijeL'fO.c), p = l,l>-\, 

x{\-\og{x))q{x)£L x {0,c), p = l,l = -i, 

for some c G (0, b). In addition, suppose H is bounded from below. If t is limit circle 
at a = we impose the usual boundary condition (corresponding to the Friedrichs 
extension; see also [4], |12j ) 

(3.37) ]im x l ((l + l)f(x)-xf'(x))=0, I £ \)- 

x — Z i 

Then 

(3.38) h = -£+ {1 + 1) V + 2) +q(x), xq & LP(0,c) 

dx z x A 



10 A. KOSTENKO, A. SAKHNOVICH, AND G. TESCHL 

and ifl>2 then 

(3.39) H =—f^ + ^—^+q(x), xqeLP{0,c). 

ax z x z 

Proof. It suffices to observe that by 28, Lem. 3.2, Cor. 3.4] 

4>(X,x) = x l+1 ui(X,x), 

where ui(X,x), xu' 1 (X,x) € W 1,p (0, c), and lim^^o xu'^x, X) = 0. Noting that 
(f>(X, x) > for all x e [0, b), by ((33|) we get 

xq(x) = xq(x) - 2x±- ^4 e iP (°' c )- 
ax ui(A, x) 

Similarly, to prove (13.39)) it suffices to note that 

9(X,x) = x~ l u 2 {X,x), 
where u 2 (X,x), xu' 2 (X,x) £ W /1,p (0,c) and \im x ^ xu' 2 (x, X) =0. □ 

Since H is limit circle at for I € [— |, h) and thus the singular Weyl function 
is a Herglotz-Nevanlinna function in this case (cf. [SHI App. A]), we obtain by 
induction: 

Corollary 3.11. Suppose 

(3.40) H = —^+ l -^^ + q(x), xq(x)eLP, 

where p G [1, oo] if I + \ tf. No and p <E (1, oo] ifl+\ G No- T/ien f/iere is a singular 
Weyl function of the form 

(3.41) M(z) = (z-A)L'+ 1 /2J Mo (z), 

where Mq(z) is a Herglotz-Nevanlinna function and X < min a(H). Here [x\ — 
max{n € Z|n < x} is the usual floor function. The corresponding spectral measure 
is given by 

(3.42) dp{t) = (t-Xi l+1 ^dp {t), 

where the measure po satisfies J R dpo (t) — oo and j R ^rff < o° ■ 

Corollary 3.12 ( 30, 28]). Assume the conditions of Corollaru \3.11\ Then there 
is a singular Weyl function from the class N%° with k = + f J • 

Proof. The inequality k < + f J follows after combining Corollary 13.111 with 
Theorem 4.2 from |50j. 

It remains to show n > \jk + |J. To this end denote by Mq(z) the Weyl function 
of the operator H with Z G [—1/2, 1/2). Then, since i? is the Friedrichs extension 
of the minimal operator in this case, Mq(z) satisfies (cf. QI3 Proposition 4]) 

(3.43) Mq(z) —> — oo, as z —> — oo. 

Moreover, by [30J Cor. A. 9] and Hypothesis 13.21 Mq{z) admits the following repre- 
sentation 

M (z) = Re(M„(i)) + / (-^ - Y^js) dpo(t), z £ [A, +oo), 
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(3.44) J 



where the measure satisfies J^°° dpoit) = oo and J^ 00 (1 + t 2 ) 1 dpo(t) < oo. Hence, 
by (|3.43p we conclude that 

r+0 ° dppjt) 

' A 1 + |*| 

Combining (|3.44|) with Corollary 13. f II and Theorem 4.2 from [30] . we arrive at the 
desired inequality. □ 

Remark 3.13. Corollary \3.12\ was first established by Fulton and Langer [15] in 
the case when the potential q(x) is analytic in a neighborhood of x = (see also 
[33] ). In the general case, it was proven in [28] (see also [30 ). Namely, by using 
high energy asymptotic of (f)(z, x) it was shown in [30] that k < [~^fcl] (for the details 
see Section 8 in [30] J . The eguality k = [5 + f J was proven in [28] with the help of 
theory of super singular perturbations and detailed analysis of solutions <p(z,x) and 
6(z,x). 

4. Connection with the double commutation method 

In this section we want to look at the effect of the double commutation method. 
We refer to [17], [21] for further background of this method. We will use the 
approach from |21j with the only difference that we include the limiting case 7 = 00 
(we omit the necessary minor modifications of the proofs of [3T], cf. [JT] Sect. 11.6]). 

Let H together with a fundamental system of solutions <p(z,x), 6(z,x) as in 
Lemma l2~2l be given. 

Hypothesis 4.1. Assume Hypothesis \2.1\ Let 7 e (0, 00] and A € R such that 
4>{X,x) satisfies the boundary condition at b if r is limit circle at b (i.e., A is an 
eigenvalue with eigenfunction <fr(\,x) if T is limit circle at b). 

Introduce 



7 + j a n\y) dy 

Denote by P and P 7 the orthogonal projections onto the one-dimensional spaces 
spanned by <f>(\,x) and (f>-y(X,x), respectively. Here we set the projection equal to 
zero if the function is not in L 2 (a,b). Note that </> 7 (A, .) G L 2 (a,b) if and only if 
7 < 00 ([U Lem. 2.1]). 

By [2TJ Lem. 2.1] the transformation 



(4.2) u 7 f(x) = f(x) - ^(A, x) / 4>{\ y)f(y)dy 

J a 

is a unitary map from (1 — P)L 2 (a,b) onto (1 — P 7 )L 2 (a, b). Moreover, by [2~T1 
Thm. 3.2] 

(4.3) H 1 {l-P 1 ) = U 1 HU- 1 {l-P 1 ), 
where the operator H 1 is associated with the potential 

(4-4) q 1 (x)=q{x)-2^\og(^- 1 + 4>{\yfdy 

and boundary conditions (if any) 

(4.5) W a $ y (A), /) = W b (0 7 (A), /) = 0. 
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Note that for 7 < 00 the operator if 7 is limit circle at a if and only H is and that 
for 7 = 00 the operator Hoc is always limit point at b. 



Theorem 4.2. Assume Hypothesis \4-l\ and let 7 < 00. The operator H 1 has an 
entire system of solutions 

<j>y(z,x) = <t>(z,x)-4>y(X,x) / (j){\ 1 y)(j){z,y)dy 

J a 

(4.6) = <j>{z, x) + — ^— -</> 7 (A, x)W x {<j>{X),cj>(z)), 

z — X 

(4.7) 9 7 {z, x) = 9(z, x) + j^-j (<£ 7 (A, x)W x {4>{X), 6{z)) + 7 7 (z, z)) . 

which satisfy W{9~ l {z) 1 <j)~ l {z)) = 1. In particular, H y satisfies again Hvvothesis \2.1[ 
Furthermore, the Weyl solutions of H 1 are given by 

(p^(z,x), ip 7 (z,x) =%/)(z,x)-\ ] —<f> 1 (X,x)W x ((f)(X),ip(z)) 

z — A 

(4.8) = 9 y (z,x) + M 1 (z)4> 1 (z,x), 
where 

Ml {z)=M{z)--\ 

Z — A 

is the singular Weyl function of H~ r The associated spectral measures are related 
via 

(4.9) dp^(x) — dp(x) +~/dQ(x — A). 

Proof. It is straightforward to check that 7 (z,x), 9~ l {z,x) is an entire system 
of solutions whose Wronskian equals one (cf. 21, (3.14) and (3.16)]). The extra 
multiple of 7 (z,a;) has been added to 9 1 {z,x) to remove the pole at z — X (cf. 
(jHH below). 

That ipry is square integrable near b follows from (cf. [2TJ (3.15)]) 

(4,o) - *(iw,<m,m)P\ 



\z-\\'dx \1 + -,S' <H\v)'d)J I 
using the Cauchy-Schwartz inequality since 

W x {cj>{X), i>{z)) = W c (cf>(X), i>(z)) + (X-z)J^ 0(A, vW(z, y)dy. 

To show that ip-y{z, x) satisfies the boundary condition at b if r 7 is limit circle at b 
we use (cf. [HJ (3.16)] plus (I3~TdD ) 

W x {cj>{X)^{z)) 



(4.11) W x (<t>y(X),^(z)) = 



i + 1 ^ a ^\v) 2 dy 

□ 



Note 

(4.12) 7 (A,.t) = 7 - 1 ^ 7 (A,a;), 7 (A, a;) = 0(A, x) + 7 (A, x)W x &(X), 0(A)). 

Remark 4.3. Again a few remarks are in order: 

(i) Clearly this procedure can be iterated and we refer to Section 4 of [21j for 
the corresponding formulas. 

(ii) If X is an eigenvalue, one could even admit 7 € [— ||(/>(A)|| -2 , 00). 
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(iii) This procedure leaves operators of the type p.35|) invariant. In particular, 
it does not change I. 

Theorem 4.4. Assume Hypothesis \4-l\ and let 7 = 00. The operator Hoo has an 
entire system of solutions 

(4.13) (j) 00 (z,x)= 1 (<p(z,x) — (ftoo(A, x) / (f>(X,y)<f>(z,y)dy 



(4.14) X)6(z,x) + <j> oo (X,x)W x (4>(X),0(z)), 

which satisfy W(9 00 (z),(f> 00 (z)) = 1. In particular, satisfies again Hypothe- 

sis[%l\ 

Furthermore, the Weyl solutions of H^ are given by 

0oo(z,x), ipoc(z,x) = (z - X)ip(z,x) + 4> 00 (X,x)W x (4>(X),ip(z)) 

(4.15) =e ca {z,x) + M 00 (z)(f> 00 (z,x), 
where 

(4.16) M 00 (z) = (z-X) 2 M(z) 

is the singular Weyl function of H^ . The associated spectral measures are related 
via 

(4.17) d Poo (t) = (t- Xfd P {t). 

Proof. In the limiting case 7 — > 00 the definition (|4.6[) from the previous theorem 
would give 0oo(A,x) = and we simply need to remove this zero. The rest follows 
as in the previous theorem. □ 

Note 

(4.18) </>oo(X,x) = 4>(z,x)-4> O0 (X,x) I 4>(X,y)4>(z,y)dy, 0^ (A, x) = (A, x) 

J a 

Remark 4.5. (i) Again this procedure can be iterated and we refer to Sec- 

tion 4 of |21j for the corresponding formulas. 
(ii) For operators of the type Q3.35P this procedure changes I to I + 2. 

5. Examples based on the generalized Backlund-Darboux 

transformation 

In this section we want to look at connections with the generalized Backlund- 
Darboux transformation (GBDT) approach (see [3SJI3I] and the references therein). 
This approach contains the double commutation method as a special case as we will 
show below and there are close relations with the binary Darboux transform (see, 
e.g., the comparative discussion in [SJ Section 7.2]). Here we will use the GBDT 
to construct an explicit example with a generalized Weyl function which is rational 
with respect to \fi and which has non-real zeros. 

More specific, we want to apply the GBDT to the Schrodinger equation 

(5.1) ru — zu, x€(0, 00). 

This case was treated in Proposition 2.2 [36], see also [35] . 

To begin with we fix an integer n £ N, two n x n matrices A, 5(0) and two 
vectors Ai(0), A 2 (0) € C" such that 

(5.2) A5(0)-5(0)A* = A(0)JA(0)*, 5(0) = 5(0)*, 
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where we have set 

(5.3) A(0) = [Ai(0) A 2 ( 



J = 



1 

-1 



Taking A(0) as initial condition we define the n x 2 matrix function A(x) 
[Ai(x) A 2 (x)] as the solution of the linear system 

(5.4) Ai(a;) = AA 2 (x) - A 2 (x)q(x), A' 2 (x) = -Ai(x) 

and S(x) via the relation 



S(x) = S(0) + / A 2 (t)A 2 (t)*dt. 



(5.5) 



Note that S(x) = S(x)* as well as the identity 

(5.6) AS(x) - S(x)A* = A(x)JA{x)* 

which follows from ()5.2|) . (|5.4|) . and ()5.5j) . Furthermore, we will assume that S(x) > 
for x > 0. 

Given these data we can construct the Darboux matrix using a transfer matrix 
function representation. To this end introduce 

(5.7) w A {z,x)=I 2 + JA(x)*S(x)- 1 (zI n -A)- 1 A{x), 

where /„ is the n x n identity matrix and the variable x is added into the transfer 
matrix function in Lev Sakhnovich form from |38) . Now one can check that (|5.7[) 
acts as a transfer matrix, that is, for any given set of linearly independent solutions 
2/0j Vi 01 (|5.1|) we obtain a set of linearly independent solutions y = [yo yi\ of a 
transformed Schrodinger equation — u" +qu = zu, where the transformed potential 
q can be expressed explicitly in terms of A and S, by virtue of 



(5.8) 



y(z,x) = [l 0] wa(z, x)w(z, x), w(z,x) 



yo(z,x) yi(z,x) 



In the special case n = 1, observe that the function A 2 is a solution of the 
Schrodinger equation corresponding to the value A of the spectral parameter. 
Therefore, in the case n > 1, the matrix A is called a generalized matrix eigen- 
value. Moreover, in the case n = 1 the GBDT contains the double commutation 
method considered in Section |4] 

5.1. The double commutation method as a special case of the GBDT. 

Let n = 1 (i.e., A is a scalar) and set 

(5.9) A = A e C, y{z,x) = [(f>(z,x) 9(z,x)], A(x) (A, x) (f>(X,x)]. 



Thus, for <p~, which is given below, formulas (|5.8I) and (|5.9p imply 



(5.10) 



<j)-r(z,x) = y(z,x) 



= [1 0] w A (z 



In view of (|5.7j) . rewrite (|5.10|) in the form 



<t>-y(z,x) = (f)(z,x) + — rr A(x) 

b(x){z — A) 



4>(z,x) 
<p'(z,x) 



4>{z, x) 
4>' {z, x) 



(5.11) 
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By f|5T2jt . (|53|) . and d5T9J) we get 

(5.12) S(x) = (A- A)- 1 A(0)JA(0)* + / \cf>{\t)\ 2 dt for A ^ A. 



In the case of A = A and real 4> treated in Section 21 identity (j5.6[) is fulfilled 
automatically. Hence, noting that 5(0) > 0, equality (|5.11[) becomes 

(5.13) <f> 7 (z, x) = <f> y (z, x) + ^^^(^(A),^)), 7 := 5(G)" 1 e (0, oo], 

2 — A 



where 
(5.14) 



7 (A,a;) = 



<f>{\x) 
S(x) 



S(x) =7 _1 + / (f>{\t) 2 dt 



Hence, if A = A and <fi(A, x) is real, the expressions for </> 7 and <f> 1 from above agree 
with gj]) and gUl, (|4J3l) . 

5.2. A generalized Weyl function with non-real zeros. For the remainder of 
this section we consider the case 



(5.15) 



0, A(0) JA(0)* = 0, 5(0) = 0. 



Clearly, (|5.2p holds for this choice of A(0) and 5(0). For q = we have (see also 

[2"4"] . where the case q = 0, 5(0) = I n was treated): 

(5.16) 

q(x) = 2((A 2 (a;)*5(x)- 1 A 2 ( a ;)) 2 + A^x)* S^)' 1 A 2 (x) + A 2 (x)* S(x)- 1 A^i)) . 



Moreover, for the case that q = 0, one can set in ([2 



(5.17) 



w(z, a;) = T{z) 



e" 1 ^ 



T(z) 



T(z) 



1 1 

u/I -\Jz 



In some places below we will need the sign in the square root fixed, hence we chose 
the branch cut in yfz along the negative real axis and assume that Im(y/z) > 0. 
Later we will need the expression for y' (see |24[ formula (1.30)]) as well: 

(5.18) y(z,x)=[~A 2 (x)*S(xy 1 A 2 {x) l] wa{z, x)w(z, x). 
If n = 1 , without loss of generality we assume that 

(5.19) A(0) = [ Vl 1] , vi = VI, 
and the second relation in (|5 . 1 5[) is immediate. 

Lemma 5.1. }30| Let n = 1 and the parameters A 6 C \ {0} and v% G R be fixed. 
Then the relation (|5.16p . where 



\A 2 {t)\ 2 dt, 



A\(x) = v\ cos(v Ae) + v Asia(V Ax), 

A 2 (x) = cos(VAx) - v\ sin(v / Ax)/\/]4, 5(x) 



explicitly determines a singular potential q satisfying q(x) = 2x 2 (1 + 0(x)) for 
x -> 0. 
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The corresponding entire solutions <p(z,x) and 6(z,x), such that both solutions 
are real-valued on R, W(6(z), 4>(z)) = 1 and <fi(z,x) is nonsingular at x = 0, are 
given by 



(z,x) = (z-A) y{z,x) 



1 

-Vi 



6(z,x) = -(z - A)y(z,x) 



The singular Weyl function corresponding to this problem has the form 
(5.20) MO) = -{z- A){z -3)(iv^+wi) _1 . 

Now, we turn to the case n = 2. To simplify calculations set 



(5.21) 



A = 



t 1 1 
fi 



, n n; A(0) = [dv v] 



d = d. 



By (|5.2ip the second relation in (|5.2[) is true and we assume also 5(0) =0. It is 
straightforward to check that (|5.4[) (with q = 0) and (|5.2ip hold for 



(5.22) Ai(x) 

(5.23) A 2 (x) 
where uj = \y/Ji, 





C3X + c 4 




C5X — C4 






Cl 




C2 


) 



(e~ uw 


C3X + c 4 




C5X — C4 




Cl 


+ e u,x 


C2 



(5.24) 



A = —\ujI2 



2u) 



1 




, A 2 = A, ci = 1 



U) 



d ci d 

C2 = 1 , C 3 = — , C 4 = — , C 5 



(5.25) 

The 2x2 matrix function S is defined by the equality 
(5.26) S(x) = f A 2 (t)A 2 (t)*dt. 



C2 

2u' 



Lemma 5.2. Let n = 2 and /ei A and A(0) 6e fixed such that (|5.2ip holds. Then 
the relations (|5.16[) and (|5.22p - (|5.26p explicitly determine a singular potential q, 
such that 



(5.27) 
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q(x) = —(l + 0(x)), x>0, x^O. 



The corresponding entire real solutions <p(z, x) and 6{z, x), such that 4>{z, x) is non- 
singular at x = and W{9{z), 4>(z)) = 1, are given by 

y(z,x) 



(5.28) (f>(z,x) 



(z - -py 



(z,x) = -(z - fi) y(z,x) 



where y is constructed via (|5.8I) . (|5.7p . and (|5.17p . Moreover, the corresponding 
singular Weyl function is given by 

(5.29) M(z) = -(z- n) 2 (z-Ji) 2 (iy/z + d)~ l . 

Proof. We begin by deriving some asymptotics at x — 0. Rewrite (|5.17p as 

(5.30) w(z, x) = cos(\/zx)l2 + sm(\/zx) ^ j- 

—\/z 
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to see that w is an entire matrix function of z (which has real- valued entries on 
and the asymptotics 

" ll 



(5.31) 



w(z,x) = (1 - -x )h +x 



-z 



+ 0{x A ) 



hold. Using (|5.22p - (l5.26l) . after some direct calculations, we get the asymptotics of 
A(x) and S(x) ±1 at x = 0. In particular, we derive 

(5.32) 

21(1 + 0(21)) x 3 {\ +0(x))' 

• 6 



T 6 AK 

detS(x) = -(l + 0(x)), 5(x)- x = - 



x 3 {\ + 0(x)) x 5 (± + 0(x)) 



45 ' x" 

Note that because of (|5.26p and the first equality in (15.32)) we have S(x) > for 
x > 0. Moreover, the asymptotics of A and formula (|5.32[) imply 

(5.33) A 2 (x)*S(x)- 1 = - [l5x- 3 (l + 0(x)) far^l + 0(x))] . 

Similarly, in view of (|5.3ip . we get 

(5.34) 

ir\ 1 x - dx 2 /2 + 0(x 3 ) x 2 /2 + 0(x 3 ) 

Hx)W[Z,X) ~[d-(z- fJ ,)x+(z-n)dx 2 /2 + 0(x s ) 1 - (z - n)x 2 /2 + 0(x 3 ) ' 

Finally, by f|5.16[) . the asymptotics of A, and f|5.32[) it is easy to show that (|5.27[) 
holds true. 

Now, consider 



(5.35) 



4>d(z,x) = y{z,x) 



1 



Clearly, <pD satisfies t</)d = z^d, where the potential is given by (|5.16j) . since 
ry = zy. According to (15.81) . (|5.7[) . (|5.21[) . and (|5.30[) . <Pd is meromorphic in z with 
only possible pole at z = /a. Moreover, we shall show that 4>d(z,x) € i 2 (0, c). For 
that purpose rewrite (|5.8[) as 

(5.36) y(z,x) = [1 0]w(z,x) + A 2 (x)*S(x)- 1 {zI 2 - Ay^Aiix) A 2 (x)]w(z,x). 



It follows from ([534]) and (zl 2 - A)^ 1 = (z _ m) 2 

1 



(5.37) {zh - A)- L A(x)w(z,x) 



z — /i 1 
z — [i 

', + 0(x 3 ) 



that 
l + 0(a; 3 ) 



d(z - fi) + 0{x) z-fi + 0(x 2 ) 



(z - 

Relations (15331 and (i05j) - (jOTj) yield <j>(z,x) e L 2 {0,c). 

To show that 4>d is an entire function, consider again the resolvent 



(5.38) {zh - A)~ 1 A(x)w(z,x) 

' A(x)w(/i, x) + 



{z-nf 



1 




A(x)w(^, x) 



(z - /i) 

d_, 

dz 



1 




A(x)w z (n,x) 





' 1 " 


) 


-d 



d 

n(z,x), 



where w z (z,x) = J^w(z,x), and fl is an entire vector function of z. In view of 
(|5.22I) - (I5.25I) and (|5.30[) . direct calculations show that the vector coefficients at 
(z — [i)~ 2 and (z — /i) _1 on the right-hand side of (|5.38l) equal zero, that is, the 
left-hand side of (15.38)) is an entire vector function. Therefore, by (|5.35l) and (|5.36[) 
one can see that <Pd(z,x) is an entire function too. 
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Let us consider another solution 

(5.39) § D (z,x) = y(z,x) !j* 

According to (|5~%|) . (|5~T71) . (|STT5|) . and ([535)1 we have 

(5.40) W(6 d (z),$d(z)) =-dctw A (z,x). 

Thus, detw,4(z,x) does not depend on x. Moreover, using (|5.7|) and (|5.6|) we get 

(5.41) uu(z,x)* Juv^z, x) = J, 

that is, I det Wa{z, x)| = 1. Taking into account the fact that it^ZjX) is a rational 
function of z with the only possible pole at z = 11 (of order no greater than 2) and 
recalling that the relations fj, ^ ~p,, 

I det wa{z, x)\ — 1, lim det wa(z, x) = 1 

hold, we derive: det wa(z, x) = {z — ~p,) k (z — [i)~ k (0 < fe < 4). Further calculations 
show that k — 2: 



(5.42) 



det wa(z, x) = (z — /i) 2 (z — fi) 



To prove that and are entire and real, rewrite (|5.41l) as Wa{z, x)* — Jwa(z, x) 1 J* 
In view of (I5.42[) the last equality yields 



w\\{z,x) u?2i(z,x) 

UJl2(z, x) W22(z, x) 



(z - /i) 2 (z - 



Wll(z,x) I0 2 l(z, x) 
Wl 2 (z,x) W 2 2(Z,X) 



where wa ='■ { w ij}i j=i- I n other words, we get 

(5.43) (z — JI)~ 2 wa{z, x) = (z — uJ^ (z, a;). 

Recall that u>(z,x) = w(z,x), that is, u> is real. Thus, by ()5.8[) and (|5.43p the 
vector function (z — /J)~ 2 y(z, x) is real for /1 7^ /J. So, the functions (/> and 9, which 
are given by ([5T281) . are real. Definitions (f5T28|) . (f5T35l) . and ([5391) yield also 



(5.44) 



(z,x) = (z-/x) 2 4>d{z,x), 9(z,x) = —{z - fi) 2 9 D (z,x). 



Since 4>d is an entire function, it follows from (|5.44[) that the real function <j> may 
have only one pole at z = ~p (p =/= ~p). Therefore, (j> is an entire function. It follows 
from (|5.44p that 9 is an entire function too. Finally, (15.401) . (|5.42p . and (|5.44j) imply 
that W(0(z), </>(z)) = 1. Thus, the statements of our lemma regarding 4> and 9 are 
proved. 

Now, using 4> a- n d we can construct explicitly a singular Weyl function M. 
Observe that (I5.17[) yields 

1 1 _ [ e [ ^ x 
is/z\ ~~ [i^/ie 1 ^ ' 

In view of (|5.23|) and (|5.26[) we calculate that 

detS'(x) ~ |c 2 c 5 | 2 (4Re(a;))- 4 e 4Rc ^) a; (x -> +00), 

and so the transfer function wa(z,x) given by (j5.7p behaves like 0(x 4 ). Hence, it 
follows from ([531) and (f5T45|) that 

1 



(5.45) 



iy(z, x) 



(5.46) ip D (z,x) = y(x,z) 



[l 6\wa(z,x) 



H/ ze 



€ Z/(c,oo). 
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Therefore, by (fOg)) . ([Q9]) . and (jQ6l) we get 

ip(z, x) — 0{z 1 x) + M(z)(j)(z, x) 6 £ 2 (c, oo), 
that is, M(z) given by (|5.29p is a singular Weyl function of our system. 



□ 



Remark 5.3. According to (|5.27|) . the GBDT generated by a 2x2 matrix A of the 
form (|5.21|) transforms a Schrodinger operator with I = into one wzt/i ^ = 3. 
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